Magnetization-based assessment of correlation energy in canted Single-Chain Magnets 
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We demonstrate numerically that for the strongly anisotropic homometallic S — 2 canted single- 
chain magnet described by the quantum antiferromagnetic Heisenberg model the correlation energy 
and exchange coupling constant can be directly estimated from the in-field-magnetization profile 
found along the properly selected crystallographic direction. In the parameter space defined by the 
spherical angles (0, 9) determining the axes orientation, four regions are identified with different 
sequences of the characteristic field-dependent magnetization profiles representing the antiferromag- 
netic, metamagnetic and weak ferromagnetic type behavior. These sequences provide a criterion for 
the applicability of the anisotropic quantum Heisenberg model to a given experimental system. Our 
analysis shows that the correlation energy decreases linearly with field and vanishes for a given value 
Her which defines a special coordinates in the metamagnetic profile relevant for the zero-field corre- 
lation energy and magnetic coupling. For the single-chain magnet formed by the strongly anisotropic 
manganese(III) acetate meso-tetraphenylporphyrin complexes coupled to the phenylphosphinate lig- 
ands, the experimental metamagnetic-type magnetization curve in the c direction yields an accurate 
estimate of the values of correlation energy A^/ks ~ 7.93 K and exchange coupling J/fcs ~ 1-20 K. 

PACS numbers: 75.50.Xx, 75.10.Jm, 75.40.Cx, 75.40.Mg 
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I. INTRODUCTION 



Molecular magnets make a wide emerging class of new 
materials whose properties extend the range of those typ- 
ically associated with magnets, i.e. they offer low den- 
sity, transparency, electrical insulation, and a possibil- 
ity of low-temperature synthesis^. Single-Chain Mag- 
nets (SCM) belonging to the class of molecular mag- 
nets, have been of particular interest since Gatteschi 
et alM- discovered slow relaxation of magnetization in 
a chain compound, comprising Co(II) centers and or- 
ganic radicals, without any evidence of phase transition 
to a three-dimensional magnetic ordering. Soon after, 
Clerac et al. discovered similar magnetic properties in 
an 5 = 3 Heisenberg ferromagnetic chain comprising 
Mn(III)-Ni(II)-Mn(III) trimers with an easy axis paral- 
lel to the chain2,. Next, other SCM systems often based 
on ferrimagnetic chains containing alternating spins of 
unequal magnitude were discovered^. The slow magnetic 
relaxation arises from large uniaxial magnetic anisotropy, 
negligible magnetic interactions between the chains and 
considerable intrachain interactionsi^. The last condi- 
tion is desirable for raising the blocking temperature. 
Moreover, for some one-dimensional quantum spin sys- 
tems with competing nearest-neighbor and next-nearest- 
neighbor interactions^ the Single-Chain Magnetic behav- 
ior can be affected with frustration^. 

The simplest description of ferromagnetic SCMs is 
based on the Glauber theory if the limit of an Ising chain 
can be explored-. The Ising ferromagnetic chains dis- 
play slow relaxation of magnetization and the relaxation 
time r depends on the spin-pair correlations. For these 
chains the correlation length ^ diverges exponentially at 



low temperatures as 



e = Coexp(A^/fcsT), 



(1) 



where is the correlation energy and is the Boltz- 
mann constant^. The exponent is proportional to the 
activation energy which enters the Arrhenius law^ and 
can be directly estimated from the experimental static 
susceptibility of the chain by the relation ^ ~ xT, plot- 
ting In(xT) as a function of 1/T. The correlation energy 
is also referred to as the energy to create a domain wall. 
Important fact is that A^ is proportional to the mag- 
netic coupling constant^ and similar conclusions remain 
valid for the antiferromagnetic Ising model^ii^ but then 
^ ^ ^/ixT). In addition, the correlation energy A^ is 
also the main part of the additive term determining the 
relaxation time in the Ising-like and anisotropic classical 
Heisenberg chains^iii. 

The relations between the correlation length and the 
zero-field susceptibility as well as between the correlation 
energy and the magnetic coupling can be used for the 
classical ferromagnetic Heisenberg model in the strong 
anisotropy limit^ii^ii^. If the Hamiltonian is defined by 
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the limit is reached when the exchange coupling J and the 
single-ion anisotropy D satisfy the relation \D/J\ > 2/3. 
Then the correlation energy calculated from the product 
xT yields the value of the coupling J because^ 

A^ = 2|J|5^ 

For ferromagnetic systems, this scenario is well 
confirmedSiiiii^ and has been exploited for a num- 
ber of compoundsi^iii, but for antiferromagnetic chains 
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with non-collinear anisotropy axeaiiii^ the situation 
is much less clear. The problem was revealed^ ^ for 
the compound of formula [Mn(TPP)02PHPh]-H20 con- 
sidered a textbook example of SCM (TPP — meso- 
tetraphenylporphyrin and PHPh — phenylphosphinate) 
and referred to as Mn-CAF (canted antiferromagnet). 
It was possible to rationalise many experimental results 
performed on Mn-CAF, assuming the generalised expres- 
sion for the correlation energy in the form 

= 2| JlS-^cosfJ, (2) 

where S is the canting angle for two classical spins ori- 
ented along the anisotropy axes. However, the propor- 
tionality between the correlation length ^ and the prod- 
uct xT was ruled out, implying that the correlation en- 
ergy cannot be determined by the magnetic measurement 
of the static susceptibility. 

In this paper we demonstrate that the correlation en- 
ergy of the quantum model of a canted antiferromagnetic 
chain can be found from the in-held single-crystal mag- 
netization profile properly chosen and its d dependence is 
given by the classical expression These results sug- 
gest that the correlation energy A^ can be estimated from 
the in-field magnetization measurements. As soon as the 
spin value, the geometrical structure and A^ are known, 
the value J of the magnetic coupling follows from Eq. 

Considering the empirical shape of the magnetiza- 
tion isotherm in the c direction for Mn-CAF, we directly 
find the J value without any fitting procedure. 

The paper is organized as follows: Sec.|lT]is devoted to 
a model describing the anisotropic quantum spin systems 
and the numerical method for accurate calculations of 
magnetic properties in a wide temperature range. Sec. lIIII 
is dedicated to presentation of results and their discus- 
sion. Sec. IIVI concludes our paper, summarizing the main 
outcomes. 



II. MODEL AND DMRG METHOD FOR 
CANTED SINGLE-CHAIN MAGNETS 

We consider the quantum anisotropic Heisenberg 
model which is needed to get quantitative estimates of 
the thermodynamic propertie o^^i^" of the chains with 
non-collinear anisotropy axes 

L-l 

i=l 

EE(^"A"'5f + A.Bi?"C''^f): (3) 

where spin S — 2 and L stands for the length of the chain. 
The exchange coupling J between nearest-neighbor spins 
is isotropic, uniform and positive. The tensors represent- 
ing the single-ion anisotropy Di and g^-factors are non- 
diagonal and depend on angles (0, 6) or {—(j), 6) for odd 



and even sites, respectively. The indices a, f3 G {a, b, c} 
define the global coordination system. The c direction is 
chosen along the chain axis and plays also the role of the 
quantization axis. 




FIG. 1. Schematic views of the chain structure in the global 
coordination system {a, b, c}, where the spherical angles {(p, 6) 
are defined. The black arrows represent the canted Mn-CAF 
spin arrangement, whereas the gray ones correspond to the 
coUinear antiferomagnetic spin arrangement {<j) = Q, 6 = 0). 

The explicit form of the tensors in the model ^ de- 
pends on the angles (</), 6) defining the anisotropy axes 
(Fig. [T]) and is given in the previous publicationsiiii^. 
We note that for all the pairs of angles {(p, 9) the canting 
angle is determined from the relation 

cos (5 = sin^(6') cos(2(/)) + cos^{e). 

Analysis of our model ([3]) is based on the numerical 
Quantum Transfer Matrix methodSSr^^ where the parti- 
tion function of the quantum chain is mapped onto the 
partition function of the classical 2d system with multi- 
spin interactions and a finite width 2M^^^. For differ- 
ent values of M, called the Trotter number, the classical 
partition functions form a series of approximants Zm, 
where the leading errors are of the order of The 
higher the value of M, the better the quantum nature 
of the problem is taken into account. As the Hilbert 
space dimension increases exponentially with increasing 
Trotter number, computations are feasible only for rel- 
atively small M. To overcome this problem, we em- 
ploy the Density-Matrix Renormalization Group tech- 
nique (DMRG)^''"'^! for determining an effective Hamil- 
tonian representation in the Hilbert subspace of a given 
size m which ensures covering the entire experimental 
temperature range. 

The DMRG method is very powerful but it is computa- 
tionally demanding as far as the resources and program- 
ming are concerned. The transfer matrices are asymmet- 
ric and the corresponding density matrices, constructed 
from the right and left eigenvectors of the largest eigen- 
value Xmax of a transfer matrix, are non-Hermitian. We 
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were able to overcome these problems by applying com- 
plex algebrai^, which guaranties the biorthogonality of 
complex eigenvectors but increases the computational 
complexity. 

In some DMRG applications^'' the Hilbert space is split 
into subspaces of fixed values of Mstag = '}2{—^ySj^ 
where Sj denotes a spin variable in a given column. 
Then the highest eigenvalue Xmax can be found in the 
block with the staggered magnetization M^tag = 0. This 
simplification follows from some symmetries which may 
occur in the four-spin DMRG vertices. For the non- 
coUinear model this is not the case, therefore \max has 
to be found in the entire Hilbert space. 

For each temperature, the free energy of the macro- 
scopic system ([3]) per site {L — >• oo) is related to the 
maximum eigenvalue \max of the transfer matrix Tlvf^, 
whereas the magnetization and magnetic susceptibility 
can be obtained by taking the first and second deriva- 
tives of the free energy with respect to external magnetic 
field. These quantities are calculated in the subspace 
Q < (j) < 7r/2 and < 6* < 7r/2 which follows from the 
Hamiltonian symmetry. 

The correlation length of the a components of spin is 
given by the following ratio^i: 



hm -^log(5f5f+«} 

it— >oo ft 



(4) 



hm hm hm - log hi'^ Im 



' M 

= hm log— — , 

M->-oo Xi, 

where S"" is the tensor product of the a-component of 
the spin operator and 2M — 1 identity matrices of 25 + 1 
size. The iV'max) vector corresponds to the Xmax eigen- 
value, whereas the eigenvalue Xi, is the highest among 
all whose eigenvector which satisfy the condition 

max 

The convergence of the approximants with respect to 
M depends on temperature and the size of the optimal 
basis set to. Our results are provided for m = 125 and 
the Trotter number M up to 20 (the lower temperature 
the higher M is needed) . According to our estimation the 
accuracy of results in the whole parameter space consid- 
ered is not lower than 0.59c3. 



III. RESULTS AND DISCUSSIONS 

Recently several compounds of antiferromagnetic one- 
dimensional systems containing Mn(III) ions have been 
reportediiii^. Here the most interesting example^^ is 
Mn-CAF, where the anisotropy axes are defined by Jahn- 
Tellcr elongation of Mn(HI) octahedra and make the an- 
gle 6 = 21.01 deg with the chain direction parallel to the c 
axis (see Fig. [IJ . Projections of the anisotropy axes onto 
the plane perpendicular to the chain, alternate along the 
a crystallographic axis making angles 4> = 56.55 deg. 



The single-crystal magnetometry measurementaii re- 
vealed three patterns of the in-field-magnetization pro- 
files. For the a direction the pattern shows linear field 
dependence and is referred to as antiferromagnetic (AF). 
The field dependencies measured along the b and c are re- 
ferred to as the weak ferromagnetic (WF) and the meta- 
magnetic type (MM), respectively. 

Taking into account the known values of (j) and 9 angles 
and performing the Monte Carlo simulations, the values 
of exchange coupling, uniaxial anisotropy and (7-factor 
were obtained (in short the Mn-CAF parameters) for the 
classical counterpart of the model ([3]) , 

J/fcs = 1.36(8) K, D/kB = -4.7(2) K, g = 1.97(1), 

(5) 

fitting the experimental susceptibility and magnetization 
curves^. So, the high value of the ratio \D/J\ ~ 3.5 
corresponding to a strong anisotropy occurred. 

Using the classical Mn-CAF parameters ([5]) and the 
corresponding angles 

= 56.55 deg, = 21.01 deg, S = 34.6 deg, 

we have determined the in-field magnetization profiles 
for three crystallographic directions. Their features agree 
qualitatively with those published previously in Fig. 5 for 
Mn-CAF—, including the sequence of patterns AF, WF 
and MM for the a, b, c directions. The MM-type profile 
for the field along the c axis was calculated earlier within 
DMRG 1^'^^ and its behavior detected some quantitative 
deviations with respect to the classical counterpart. 
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FIG. 2. The angle-dependent diagram showing sectors with 
different sequences of the in-field magnetization patterns. In 
sectors A, B, C the MM-type patterns occur in the a, b, c 
directions, respectively. All the patterns are outlined in the 
inset in the zero-temperature limit. 

Relaxing the angles 4> and 9, we have noticed that the 
sequence of the magnetization patterns varies with the 
easy axes orientation and we have established the do- 
mains (Fig. in the (0, 9) parameter space. Within 
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a given domain the fixed shape of the low-temperature 
magnetization isotherms as a function of field for all the 
crystallographic directions occurs. The expected zero- 
temperature patterns are plotted in the inset. The sec- 
tors appearing in the diagram are described in Table ID 
In sectors A, B, C the MM-type patterns are unveiled 
for the field along the a, 6, c axes, respectively. The 
sequence of the AF and WF profiles accompanying the 
MM-type shape in a given sector is specified in Table HI 
In the B sectors the AF pattern is missing and two WF 
patterns are present in return. The crystallographic di- 
rection, where the magnetization values are higher in the 
low fields, has been distinguished by an additional index 
a (or 6). 

TABLE I: Pattern sequences in the diagram sectors. 



Sector 


Crystallographic 


direction 


symbol 


a 


b 


c 


A 


MM 


WF 


AF 


B, 


WF 


MM 


WF 


Be 


WF 


MM 


WF 


C 


AF 


WF 


MM 



The B sectors in the diagram are separated from the 
other sectors by the dashed line obtained from the so- 
lution of the equation cot(26') — sin(2(/)) cos(20). This 
line corresponds to the canted angle between the adja- 
cent anisotropy axes which amounts to S = tt/2. The 
solid line in the diagram follows from the symmetry of 
the Hamiltonian and can be determined from the solution 
of the equation cot 9 = cos 0. 

The diagram in Fig. [2] provides the criterion whether 
the model ^ can describe a particular magnetic system. 
It predicts the occurrence of the MM-type profile which 
is accompanied by a given combination of the AF-type 
and the WF-type patterns. 

We immediately see that the criterion is fulfilled for 
Mn-CAF— . Its geometric coordinates (cj), 9) belong to 
the sector C and the experimental single-crystal profiles 
agree with the patterns predicted for this sector. 

As the zero-field susceptibility and the correlation en- 
ergy are not related for the canted antiferromagnetic 
chains, we have decided to analyse the field dependence of 
the correlation energy (H). First we checked if the cor- 
relation length diverges in the limit of low temperatures. 
For low temperatures the upper panel of Fig. [3] demon- 
strates a linear dependence of In(^c) (where is deter- 
mined from Eq. ([4])) on the inverse of temperature, choos- 
ing the Mn-CAF parameters ^ and the corresponding 
angles = 56.55 and 6 = 21.01. This dependence implies 
that the low-temperature correlation length ^ exhibits 
the exponential divergence in the MAf -direction 

C^CoeMMH)/kBT), (6) 

but the slope A^(iJ) is field dependent. 
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FIG. 3. (Above) Temperature behavior of the correlation 
length for various values of the magnetic field H parallel to 
the MM axis (the Mn-CAF parameters). (Bellow) Our esti- 
mations of the A({H) vs external magnetic field parallel to 
the c axis plotted by triangles. The solid line represents the 
linear regression: /(_ff) = 9.10 — 2.33 * -ff. The magnetization 
isotherms at T = 1 and 1.6 K are given by the dashed and 
chain lines, respectively. 

The field dependence of the correlation energy A^{H) 
has been studied numerically and plotted by triangles in 
the bottom part of Fig ([3|) . We have revealed the striking 
linear behavior which can be very well reproduced by the 
function 

A^{H)/kB = 9.10-2.35* H, (7) 

where the right hand side is expressed in K and the field 
is given in T. So, the function Aj(iJ) vanishes, which im- 
plies a finite value of correlation length ^ (see Eq. ([S])), for 
the critical field Her = 3.91 T, what is close to the value 
estimated directly from the DMRG correlation length. 
Her = 3.86 T. 

The abrupt change of the correlation length from the 
infinite to the finite value induced by the field should 
be accompanied by some changes in the physical prop- 
erties. We demonstrate in the bottom panel of Fig. [3] 
that the MM-type profile displays a rapid jump in the 
vicinity of Her T from to a nearly saturation value. We 
interpret this behavior as a field-driven transition from 
the antiferromagnetic to the ferromagnetic order along 
the c direction. Such critical field can be determined 
from the intersection of the low-temperature magnetiza- 
tion isotherms, from their infiection point or from the 
maximum of the in-field magnetic susceptibility. 

Now we analyze the correlation energy Aj for the Mn- 
CAF parameters in the absence of a magnetic field. The 
formula implies A^{H = 0)/kB = 9.10 K, while the 
DMRG calculations based on Eqs. ([T]) and ^ yield 
A^/ks = 9.05 K. These estimates are in remarkable 
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agreement with the classical prediction^ ^ = 8.95 K 
following from Eq. ([5]). Moreover, the coincidence of 
the correlation energy of the quantum model ([3]) and the 
classical dependence ([2]) is confirmed and demonstrated 
for a number of couplings J and the canting angles 5 in 
Fig.H 
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FIG. 4. Comparison between the functions 25^1 cos((5)| and 
5''^|cos(5)| vs 5 plotted by the full and the dashed line, re- 
spectively, and the ratios A^/J and ^bMctHct / J illustrated 
by the circles and squares. The open symbols correspond to 
J/fcs = 1.26 K and the full symbols to J/ks ~ 1.36 K. 

In Fig. 13] the upper continuous and lower dashed lines 
represent the functions 2S^ cos 6 and cos 6, respec- 
tively. The open and full circles illustrate the ratios A^/J 
calculated numerically within DMRG for J/fcB = 1.26 K 
and J/A:b=1.36 K, respectively. As to the angles 6, their 
values follow from the angles (j) and 9 chosen from dif- 
ferent sectors of the diagram in Fig. [21 The fields were 
applied in the proper direction for a sector selected. 

We conclude that the agreement between the classical 
predictions 2S^cos(5 for the ratio A^/J and those calcu- 
lated for our quantum model ^ and plotted in Fig. 2] is 
very good. There are some deviation for small and high 
values of S, but they can be attributed to uncertainties of 
the extrapolations performed for In(^) yielding the values 
A J . These uncertainties are of the order of a few percent 
and are much higher than the accuracy of the DMRG 
results. They could be diminished lowering temperature 
and increasing the cost of our simulations. Concluding, 
our numerical results in Fig. 2] give a strong evidence 
that the correlation energy of the canted antiferromag- 
netic quantum chains can be expressed by the relation 
([2]) previously found for the classical systems. 

In Fig. [4] we also plot by the full squares the ratios 
fisMcrHcr/ J as a function of S, where Mcr is the value 
of the MM-type magnetization profile corresponding to 
the field Her- We emphasise that the ratio is com- 
pletely determined by the coordinates of the inflection 
point of the MM-type magnetization profile. The ratios 
fisMcrHcr/ J calculated coincide with the values of the 
function cos 6 drawn by the dashed line and are lower 



by a factor of 2 than the ratios A^/J, implying that 

Aj = 2nBMcrHcr- (8) 

We can interpret the result ([5]) that the correlation energy 
is equal to the average Zeeman energy per pair of spins 
calculated for the critical field Her- 

The relation ([5]) is an important result of our paper 
because it provides a direct estimate of the correlation 
energy Aj in terms of the special coordinates of the 
MM-type magnetization profile. Moreover, the exchange 
constant J/fc^ can be also obtained from the MM-type 
single-crystal low-temperature magnetization isotherm. 
Combining Eqs. and ([5]), we find the formula for the 
magnetic coupling 

J = flsMcrHcr/iS'^COsS). (9) 

To validate our conclusions, in Fig. [S] the MM-type 
behavior is drawn for the ratios \D/ J\ ~ 3.5 and \D/ J\ ~ 
4.5 in both panels. For the angles specified in the legends 
and coming from the sectors A and B, the a and b panels 
display the profiles found for the magnetic field parallel 
to the a axis and b axis, respectively. The slope of all 
the MM-type curves in Fig. [S] increases with decreasing 
temperature and the shape does not change. 

The coordinates of the intersections of the correspond- 
ing magnetization isotherms in Fig.[5](T = 1.0 K and 1.6 
K) inserted into Eq. ^ yield the consistent estimates of 
J/ks with the uncertainty equal to ±0.04 K for the input 
values J/ks — 1.26 K and 1.36 K. When we used lower 
temperatures (T = 0.5 K and 1.0 K), the uncertainty 
declined to ±0.005 K for the same J/fc^, which is the 
exDCcted trend. 




1 2 3 4 5 6 

H[T] 



FIG. 5. Metamagnetic profiles as a function of magnetic field 
for two sets of parameters presented in the legends. The thick 
lines correspond to T = 1.6 K and the thin lines to T = 1.0 K. 
The solid and the dashed lines belong to the A sector {(p — 10 
deg, — 80 deg) and to the Be sector {(j) — 75 deg, ^ = 65 
deg), respectively. 

Now, taking the values S = 34.6 deg, H^r = 3.45 T 
and Mcr = 1-71 (in (Ib) directly from the experiment for 
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Mn-CAF compoundii, we can evaluate from Eqs. ([8]) and 
© the value of correlation energy A^/ks = 7.93 K and 
the coupling constant J/fc^ = 1.20 K, respectively. Both 
the values are slightly lower than the A^/ks — 8.95 K 
and J/kB = 1.36(8) K reported earlier^^. However, this 
fact can be understood as the experimental value J/ks = 
1.36(8) K was settled on the basis of the classical model 
and the fitting procedure which usually overestimate the 
coupling constant^. If the coupling J is overestimated 
then the relation Q may proliferate the inaccuracy and 
imply somewhat higher A^. 

IV. CONCLUSIONS 

Our investigation has confirmed that the relative po- 
sitions of the anisotropy axes affect the shape of mag- 
netization curves along different crystallographic direc- 
tions. The resulting angle-dependent diagram predicts 
four areas, each characterized by a specific sequence of 
the in-field magnetization patterns. If for a compound 
studied the shapes of magnetization curves in the low- 
temperature range are consistent with those established 
in the diagram for a given geometry, then the compound 
can be analyzed in terms of the quantum Heisenberg 
model with a strong anisotropy. 

For the canted antiferromagnet SCM the values of cor- 
relation energy (H) decrease linearly with the applied 
field and vanish at Her which coincides with the intersec- 
tion point of the MAf-type magnetization profiles. The 
coordinates Her and A/cr of this special point determine 
the zero-field correlation energy A^ which in turn reveals 



the classical dependence on the canting angle S and the 
coupling J. These findings relate the J value to the co- 
ordinates Her and Mcr and yield a simple recipe to the 
evaluation of J. 

Following these ideas, we have estimated for the Mn- 
CAF compound the relevant quantities at A^/ks — 7.93 
K and J/ks = 1.20 K which agree very well with the 
values reported earlier. 

Our calculations provide some insight into the common 
claim that in the strong anisotropy limit the Heisenberg 
model can be substituted by the Ising model. We confirm 
that the projections of spins onto the axis displaying the 
MM-type pattern can be treated as classical Ising-type 
variables. Then the rapid jump in the magnetization 
profile can be related to formation of an Ising domain 
wall. Anyway, it is striking that the quantity Aj defined 
in the absence of field can be so easily settled from the 
experiment performed in the presence of field. 
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